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ABSTRACT 

We report the discovery of low-amplitude gravity-mode oscillations in the massive binary 
star V380Cyg, from 180d of Kepler custom-aperture space photometry and 5 months of 
high-resolution high signal-to-noise spectroscopy. The new data are of unprecedented quality 
and allowed to improve the orbital and fundamental parameters for this binary. The orbital 
solution was subtracted from the photometric data and led to the detection of periodic intrinsic 
variability with frequencies of which some are multiples of the orbital frequency and others 
are not. Spectral disentangling allowed the detection of line-profile variability in the primary. 
With our discovery of intrinsic variability interpreted as gravity mode oscillations, V380Cyg 
becomes an important laboratory for future seismic tuning of the near-core physics in massive 
B-type stars. 

Key words: binaries: eclipsing — stars: individual (V380Cyg) — stars: fundamental param- 
eters — stars: variables: general — stars: oscillations 



1 INTRODUCTION 

V380Cyg (HD 187879, KIC5385723) is a bright (V = 5.68) 
double-lined detached eclipsing binary with two early-B type com- 
ponents in a 12.426 d eccentric (e = 0.23) orbit. It is an important 
test for massive star models, with a primary star about to start the 



* Based on DDT data gathered with NASA's Discovery mission Kepler and 
with the HERMES spectrograph, installed at the Mercator Telescope, oper- 
ated on the island of La Palma by the Flemish Community, at the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astrofi'sica 
de Canarias and supported by the Fund for Scientific Research of Flan- 
ders (FWO), Belgium, the Research Council of K.U. Leuven, Belgium, the 
Fonds National de la Recherche Scientific (F.R.S.-FNRS), Belgium, the 
Royal Observatory of Belgium, the Observatoire de Geneve, Switzerland 
and the Thiiringer Landessternwarte Tautenburg, Germany. 



hydrogen shell burning phase and a secondary still in the early part 
of the main seque nce. V380 Cyg has been t h e subject of numerou s 
stu dies, of which IPopper & Guinart dl998l) . iGuinan et all | |2000|) . 
and lPavlovski et alj J2009h are the most extensive, to which we re- 
fer for details on the properties of the system. 

The inteipretation of the stellar parameters of V380 Cyg de- 
duced from the binary data is problematic in the sense that an ex- 
treme core overshooting parameter value above 0.5 (in local pres- 
sure scale height s) is needed to match them to theoretical evo- 
lutionary tracks dPavlovski et al.ll2009l) . Such a high core over- 
shooting is not supported by recent seismic estimates, which rather 
point to values below 0.2 for single B- t ype pulsators with simi- 
lar masses (e.g.. lAerts et al.ll2003l 1201 it lPamvatnvkh et alj|2004 



lar masses (e.g.. lAerts et al.ll2003L 1201 U lPamvatnvkh et al.1120041: 
iBriguet et alj201 lh . Nevertheless. lBriquet et alj | |2007l) deduced an 
overshooting parameter of 0.44 ± 0.07 for the 8.2 ± 0.3 M Q primary 
of the detached spectroscopic binary 6 Oph. 



2 A. Tkachenko et al. 



.1 

-A 

\ 

■i 


1 

1 1 

1 


1 1 

\ A 1 

111 

I 

1 1 


1 

i 

\(\ 

\ 

1 


1 

\a! 

! 

i 


i 

f 

i 


1 

ill 

l 


1 1. 

1M 

- 


( if w . 


i i 

v ' v V 


H b: 

< ; 

% \ 


1 1 1 1 
■ 1 1 1 1 


1 1 
1 1 


i i i i 

■ t ' i ■ 

i i i i- 



10 20 30 40 50 60 70 80 

BJD - 2455641.510101 



Figure 1. Upper panel: Kepler light curve of V380Cyg observed during Q9 (gray) with the binary model overplotted (black). Middle panel: residual light 
curve. Lower panel: first moment of the disentangled Si III line of the primary. 



Clearly, an independent seismic evaluation of the core over- 
shooting parameter of the components of V380 Cyg is of major im- 
portance. However, despite the position of the two binary compo- 
nents in the instability strips for B stars, searches for oscillations in 
the binary so far led to null detections. We report the discovery of 
intrinsic variability from uninterrupted high-precision space pho- 
tometry obtained by the Kepler satellite and interpret it in terms of 
gravity (g-mode) oscillations of the primary. 



2 OBSERVATIONS 

2.1 Kepler customized aperture photometry 

The case study of the prototypical pulsator RRLyr, with a maxi- 
mum brightness of V = 7.06 during its pulsation cycle, showed that 
Kepler's instrument is capabl e of performing high-precision pho- 
tomet ry on saturated targets dBrvson et al.l boiO; Kol enberg et al .1 
l20l ]]) by using a customized aperture mask. A dedicated mask was 
thus defined to collect all flux from V380 Cyg and the binary was 
successfully observed in short cadence (58 s) mode during Quarter 
7 (Q7) of the mission, for about 90 d continuously. This led to the 
conclusion that low-amplitude variability occurs, aside variability 
from the eclipses and brightenings, so the star was reobserved dur- 
ing Quarter 9 (Q9). 

The Q9 Kepler light curve is shown in Fig.[T] Variability on 
timescales shorter than the orbital period is obvious throughout the 
entire dataset, including during the eclipses. The Kepler light curve 
is fully compatible with but more precise than ground-based pho- 
tomet ric data, which had previously revealed scatter at the mmag 
level l lGuinan et al.l2000l) . 



2.2 HERMES high-resolution spectroscopy 



Given the discovery of low-amplitude variability in the Kepler Q7 
data, we initiated a ground-based spectroscopic campaign. We ac- 
quired an extensive time-series of high-resolution, high signal-to- 
noise (S/N) spectra, starting during Kepler Q9 and lasting for 142 
days. The spectra were taken with HERMES, the fibre-fed high- 
resolution spectrograph on the Mercator telescope (Observatorio 
del Roque de los Muchachos, La Palma, Canary Islands), which 
samples t he entire optical ra nge (380-900 nm) with a resolution 
of 85 000 dRaskin etal]|201il) . We obtained 150 observations dur- 
ing out-of-eclipse phases and 256 exposures during primary or sec- 
ondary eclipses. The typical exposure time was 1200 s. The data 
taken simultaneously with the Kepler Q9 observations is repre- 
sented on Fig.Q] 



The reduction of the spectra has been carried out using the 
dedicated HERMES software pipeline, including bias subtraction, 
cosmic ray filtering, flat fielding, wavelength calibration using a 
ThAr lamp, and order merging. Normalization to the local contin- 
uum was done manually by fitting a spline function to carefully 
selected continuum points. We c omputed Least Squa res Decon- 
volution (LSD) profiles following iDonati et al .1 dl997h . These are 
shown in a gray-scale re presentation accord ing to the orbital period 
of P orb = 12.425719 d jGuinan et al1l2000|) in Fig.|2] The double- 
lined nature of the binary is clearly recovered in our high-quality 
spectroscopy. 
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Table 1. Summary of the parameters for the wd2004 solution of the Kepler 
short-cadence light curve of V380 Cyg. For f urther details on the control pa - 
rameters we refer to the wd2004 user guide (Wilson & Van Hamme 2004). 
A and B refer to the primary and secondary stars, respectively. 



300 



Parameter 



wd2004 name Value(±Error) 



Control parameters: 

wd2004 operation mode mode 

Treatment of reflection mref 

Number of reflections nref 

Limb darkening law ld 

Numerical grid size (normal) n1,n2 

Numerical grid size (coarse) nIl, n2l 

Fixed parameters: 

Orbital period (P 01 -b, d) period 

Reference time of minimum (BJD) hjdO 

Mass ratio (q) rm 

Effective temperature A (T^, K) tavh, tavc 

Effective temperature B (r?, K) tavh, tavc 

Rotation rates f1, f2 

Gravity darkening gr1,gr2 

Bolometric albedos alb 1 , alb2 

Bolometric LD coeff. A xbolI 

Bolometric LD coeff. B xbol2 

Third light el3 

Passband LD coeff. B x2 

Fitted parameters: 

Phase shift pshift 

Star A potential phsv 

Star B potential phsv 

Orbital inclination (/, °) xtncl 

Orbital eccentricity (e) e 

Longitude of periastron (a), ° ) perrO 

Light from star A hlum 

Light from star B clum 

Passband LD coeff. A x2 

Derived parameters: 

Fractional radius of star A (R^/a) 

Fractional radius of star B (Rq /a) 

Physical properties: 
Mass of star A (Ma, M g ) 
Mass of starB (Ma, M q ) 
Radius of star A (R A , R G ) 
Radius of starB (R B , #o) 
Surface gravity star A (logg) 
Surface gravity star B (logg) 
Orbital semimajor axis (a, R e ) 

log(L A /L ) 
log(L B /L ) 
Distance (d, pc) 





2 (detailed) 
2 

1 (linear) 
60 
40 

12.425719 
2441256.544 
0.62 
21500 
22000 

1.0 

1.0 

1.0 
0.648 
0.685 

0.0 
0.262 

-0.0493±0.0003 
4.62±0.05 
11.1 ±0.6 
80.5±0.5 
0.235±0.005 

135+2 
12.02±0.10 
0.777±0.011 
0.08±0.03 

0.2655±0.0020 
0.0648±0.0011 

11.80±0.13 
7.194±0.055 

16.00±0.13 
3.904±0.067 
3.102±0.007 
4.112±0.015 

60.25±0.19 
4.711±0.023 
3.474±0.047 
987±17 



3 IMPROVED ANALYSIS OF THE BINARITY 
3.1 Spectral disentangling 

We applied the technique of spectral disentangling (spd) in velocity 
(Four ier) space as implemented in the FDBbinary code dlliiic et all 
120041) . using three spectral intervals centred on strong helium or 
metal lines (He I 4471 A and Mgll 4481 A; He I 4920 A; the Si III 
triplet at 455CM1575 A). The orbital elements we deduced are the 
time of periastron, r pcri , the eccentricity, e, the longitude of peri- 
astron, a), and the amplitudes of the RV variations for the primary 
(star A) and secondary (starB) components, A" A and Kb. 

Since the effective temperatures (T c g) of the two stars are sim- 
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Figure 2. Gray-scale representation of the LSD profiles folded according 
to the orbital period of 12.4257 19 d. Phase ze r o corr esponds to the time of 
primary minimum as derived bv lGuinan et al 1 d2000l) . 



ilar dGuinan et al 1 l2000l : I Pavlovski ct al ] l2009h . their spectra are 
also similar at optical wavelengths. The orbital solutions for all 
three considered short spectral intervals led to a consistent set 
of values: r pcn = 54602.824 ± 0.161 d, e = 0.2338 + 0.0026, 
oj = 139.54 ± 0.58 degrees, K A = 94.266 ± 0.034 km s" 1 and 
A"b = 154.59 + 0.31 km s . The quoted errors are the standard devi- 
ations of the spd solutions for the three selected wavelength ranges. 



3.2 Light curve modelling 

The orbital effects in t he light curves were mod elled using the 
Wilson-Devinney code dWilson & Devinnevlll97ll) in its 2004 ver- 
sion (hereafter wd2 004), with automatic ite ration performed using 
the jktwd wrapper dSouthworth et al.l201 ll). We converted the data 



to orbital phase using the ephemeris of lGuinan et alj ( 1200 0), sorted 
them and reduced them into 202 binned datapoints. The bins cover 
the eclipses five times finer than other phases, but are otherwise 
equally spaced. 

A variety of possible solutions were tried, with different num- 
bers of fitted parameters and model options. The full set of input, 
output and control quantities for our final solution is given in Ta- 
bleland the fit is shown as a full line in the upper panel of Fig.[T] 
The fit was performed in Mode 0, which renders the effective tem- 
peratures inconsequential by decoupling them from the solution. 
The detailed treatment of reflection was found to be necessary for 
V380 Cyg. We assumed gravity darkening coefficients and albedos 
of unity. Third light was found to be negligible so was set to zero. 
A linear limb darkening law was used, and bolometric coefficients 
for both stars plus the AT p -band co efficient for starB w ere fixed to 
values obtained using the tables of IVan Hammej Jl993h . The uncer- 
tainties on the photometric parameters listed in TableQ]come from 
the scatter of appropriate solutions with different sets of fitted pa- 
rameters, and are much greater than the formal errors of the best 
single solution. 

For our adopted solution we fitted for a phase shift with re- 
spect to the orbital ephemeris, the potentials and X" p -band light 
contributions of the two stars, the A" p -band limb darkening coef- 
ficient of the primary star (star A), e, a> and the orbital inclination. 
These parameters are given in the lower part of Table[T]alongside 
uncertainty values obtained by comparing all preliminary solutions 
providing a good fit with physically reasonable parameter values. 
We also give the fractional radii (stellar radii divided by the orbital 
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Time [d] Amplitude [mmag] 

Figure 3. Left: Short Time Fourier Transform of the residual light curve shown in the middle panel of Fig.[T] for a window width of P or b- Right: summed 
amplitude spectrum. Thick white lines near the top of panels indicate the position of periastron passages. 



semimajor axis) which are needed to calculate the physical proper- 
ties of the tw o stars. There is a goo d agreement with the parame- 
ters found bv lPavlovski et al ] ll2009h . Table[Qalso contains the full 
physical pro perties of the V380Cyg system, calculated using the 
absdim code l lSouthworth et al. I l2005h . our new sp ectroscopic and 
photom etric results, and the T^s values found by IPavlovski et al] 
J2009h . For distance we quote the value obtained usin g the 2MASS 
K s app arent magnitude, bolometric corrections from bessell et al.l 
dl998h . and adopting a reddening of E(B - V) = 0.21 ± 0.03 mag. 

We estimated the projected rotational velocity of V380 CygA 
from a spectrum taken at periastron and one at apastron, taking into 
account rotational and thermal broadening but ignoring pulsational 
broadening — see Sect.|4] The values are 91.5±7.0kms 1 and 
96.5+7.0 km s -1 , respectively, and point to a higher rotation rate at 
apastron than at periastron, although the two values are within their 
error bars. In any case, none of the two values is compatible with 
synchronous rotation, which would require 66. 1 ± 1 . 1 km s for the 
radius we found. Assuming the rotational axis to be perpendicular 
to the orbital plane, not necessarily the case in an eccentric non- 
synchronised binary, we estimate f mt to be 0.1 15 d~' at periastron 
and 0.121 d~' at apastron. 



4 THE INTRINSIC VARIABILITY OF THE PRIMARY 

The residuals of the light curve of both Q7 and Q9 (see middle 
panel of Fig.Q] for Q9) were subjected to frequency analysis af- 
ter detrending. The outcome in the region [0, 4] d~' ([0, 46] yuHz) is 
shown in the right panel of Fig.[3] Strong excess power is present 
at numerous frequencies below 1 d~' (1 1.57 /jHz), the highest am- 
plitude signals being very significant — see Table|4] where the sig- 
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wavelength (Angstroem) 

Figure 4. Snapshots of the Si III line and its disentangled version for star A 
shifted to its center of gravity after subtraction of the disentangled profile 
of the secondary. 



nificance level was computed from the local average amplitude in 
the spectrum over a range of 3 d -1 . In addition, isolated peaks oc- 
cur in [1, 4] d - '. Most of the significant frequencies are multiples of 
the orbital frequency, while others are not. The Short Term Fourier 
Transform computed for a window width taken equal to the orbital 
period (left panel of Fig.[3j shows that part of the variability damps 
and re-appears throughout the data set. 

We also searched for residual variability in the strong Si III line 
of the primary after spd (Fig.[4]l. We computed the line moments 
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Table 2. Results of the frequency analysis from the various diagnostics af- 
ter orbital subtraction. The Rayleigh limit amounts to 0.0036 d _1 for the 
photometry and 0.0071 d - ' for the spectroscopy. 



Diagnostic 


Frequency 


Factor 


Amplitude 


Significance 




td -1 ) 


forb 




level (S/N) 


Kepler LC 


0.2572 


3.20 


1.40 


8.0 


(mmag) 


0.0802 


1.00 


1.25 


7.1 




0.1295 


1.61 


1.15 


6.6 




0.2769 


3.44 


1.02 


5.8 




0.3204 


3.98 


1.00 


5.7 




0.1611 


2.00 


0.98 


5.6 




3.4611 


43.01 


0.16 


3.6 




3.3802 


42.00 


0.16 


3.5 




2.2535 


28.00 


0.31 


3.5 




0.8053 


10.01 


0.56 


3.2 


< v > 


0.0800 


0.99 


1.64 


8.5 


(kms~') 


0.8062 


10.02 


1.17 


6.1 




0.8853 


11.00 


0.94 


4.9 




0.2738 


3.40 


0.80 


4.2 




1.3682 


17.00 


0.57 


3.0 


Pixel-by-pixel 


0.8868 


11.02 


0.0104 


5.2 


(continuum units) 


0.0350 


0.43 


0.0118 


3.8 




0.0797 


0.99 


0.0099 


3.2 



following the definition of lAerts et al. I d2010l) of all 406 spectra, 
corrected for the disentangled profile of the secondary and shifted 
to the center of gravity of the primary. The (v)-values obtained in 
this way for the spectra taken simultaneously with the Kepler Q9 
photometry are shown in the lower panel of Fig. Q] while the fre- 
quencies found in this indepen dent data set are also listed in Ta- 
ble|4] A pixel-by-pixel analysis (;Zima 2008) was also done. Com- 
paring all frequency analyses results in Table|4] leads to the con- 
clusion that the line-profile variability of the primary is compat- 
ible with the Kepler residual photometry in the sense that power 
is found at multiples of / olb and at some other frequencies. The 
frequencies in the disentangled Si III lines cannot be due to the 
prewhitening of the orbital motion, because the variability occurs 
as asymmetries in the profiles and are thus intrinsic to the primary 
(FigS. 



5 DISCUSSION AND INTERPRETATION 

The intrinsic photometric and line-profile variability detected in the 
residual data of V3 80 Cyg A after orbital subtraction are naturally 
explained in terms of gravity-mode oscillations. Other causes of 
the intrinsic variability, such as migrating spots or rotational mod- 
ulation, are unlikely as they would require a complex time- varying 
surface pattern and one would expect a much more prominent pres- 
ence of the rotational frequency in that case. We also found no evi- 
dence of mass loss in the spectroscopy. 

Some of the intrinsic frequencies detected in V380 Cyg A are 
connected with the orbital frequency. Given their time-dependent 
amplitudes (Fig.[3j we suspect that they are tidally triggered at pe- 
riastron while damped further in the orbit. If confirmed, this finding 
holds the potential to model the interior structure of the primary in 
great detail and to deduce a seismic estimate of the core overshoot- 
ing parameter which can then be compared with the one deduced 



from isochrone fitting. However, additional data are needed to de- 
duce if the amplitude and frequency behaviour is indeed repeti- 
tive for numerous orbits. Such a situation was found for few close 
binaries so far. One case is HD 174884, an eccentric short-period 
binary consisting of two cool B stars with residual variability at 
low-order multiples of the orbit al frequency and two tidal oscil- 
lations induced at 8 and 13 f m h (M aceroni et al.|[2009h - Resonant 
mode locking at 90 and 9 1 / or b occurs in the highly e ccentric 42 d 
binar y HD 187091 consisting of two A stars (KOI-54. IWelsh et al.l 
1201 ll) . In contrast to these two cases, V380 CygA exhibits an excess 
bump of power rather than only well-isolated stable frequencies. 
Future Kepler data will cover three more quarters while additional 
spectroscopy will be gathered to unravel the pulsational behaviour 
of this massive binary. 
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